Observation of a two-dimensional Fermi gas of atoms 
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We have prepared a degenerate gas of fermionic atoms which move in two dimensions while 
the motion in the third dimension is "frozen" by tight confinement and low temperature. In situ 
imaging provides direct measurement of the density profile and temperature. The gas is confined in 
a defect-free optical potential, and the interactions are widely tunable by means of a Fano-Feshbach 
resonance. This system can be a starting point for exploration of 2D Fermi physics and critical 
phenomena in a pure, controllable environment. 



Two-dimensional Fermi systems are predicted to have 
rich physics of phase transitions and quantum critical 
points |l|, |2J ■ Reduction of the spatial dimensionality in- 
creases the role of fluctuations, which in turn causes phe- 
nomena such as superfluidity without Bose-Einstein con- 
densation and non-Fermi-liquid behavior. Reduced di- 
mensionality along with strong interactions is also at the 
origin of high-temperature superconductivity. 2D collec- 
tive phenomena have been studied in helium-3 films |3[ 
and in a rich collection of systems containing electron gas 
in superconducting [J, |5| and nonsuperconducting [y, |7| 
phases. 

An ultracold atomic Fermi gas may provide unique 
possibilities for studying 2D phenomena in a controlled, 
impurity- free environment. In 3D atomic Fermi gases |8|, 
it has been possible to tune at will the principal exper- 
imental parameters such as interactions [9(, energy [lOj, 
and spin composition llj, [12( and to have more than 
two spin states [13] . This control and tunability un- 
precedented to other Fermi systems should also be avail- 
able in perspective experiments with the 2D ultracold 
gases. Experiments with 3D optically trapped atomic 
Fermi gases 1J, LL5| and their derivatives have provided 
first-time observation of fundamental quantum phenom- 
ena, which include superfluidity 16l4l9J| and mechanical 
stability [20( of a resonantly-interacting Fermi gas, co- 
herent transformation of a Fermi system to a Bose sys - 



21 23], 



tern |9[ , Bose-Einstein condensation of molecules 

and possibly the shear viscosity near its fundamental 

quantum minimum 2J, [25[. Three-dimensional ultra- 



cold Fermi gases have also been used to test theoretical 
models of other Fermi systems: neutron stars and nu- 
clear matter 2fJ, |2jj , quark-gluon plasma [28( , and high- 
temperature superconductors [29|. Availability of a 2D 
atomic Fermi gas would allow one to study Fermi mix- 
tures in mixed dimensions, where one species is confined 
to 2D while the other to 3D 0, [Il|. In the field of 
ultracold gases, Bose systems are experimentally avail- 
able in three, two, and one spatial dimension [32] . The 
Fermi systems, however, have been created only in 3D 
and ID [8|, |32|, while 2D degenerate Fermi systems have 
been missing. 

In this Letter, we report on preparation and in situ 



imaging of a 2D atomic Fermi gas. In situ observation 
of a density distribution is among the most direct meth- 
ods for probing a quantum system. This capability of 
an atomic gas is unique to presently available 2D Fermi 
systems. It may let one see phase separation as well 
as measure thermodynamic, statistical, and mechanical 
properties, the amount of the mean field, and in some 
cases, the phase of a many-body wave function. 

Two-dimensional kinematics is obtained by confining 
atoms in a highly anisotropic pancake-shaped potential 
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and keeping their energy below the energy of the first 
axial excited state. 

A series of pancake-shaped potentials is created by set- 
ting up a standing optical wave along the z direction as 
shown in Fig. [TJ The frequency of the standing wave is 
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FIG. 1. Trapping ultracold atoms in antinodes of a standing 
optical wave. The isolated clouds of atoms shown in dark red, 
the standing-wave intensity shown in light purple. 

far below the resonance of the atoms. As a result, the 
minima of the dipole potential are at the intensity max- 
ima and each antinode acts on the atoms nearly as the 
potential of Eq. (JT]). The Gaussian shape of the mode as- 
sures weak confinement along the transverse directions. 
In Fig. [2] one may see a snap shot of the density distri- 
bution taken along the y direction, parallel to the plane 
of the pancake-shaped clouds. The gas is prepared in 
a strongly degenerate regime, with the temperature of 




FIG. 2. In situ image of the column density distribution 
ni{x,z) of the Fermi gas trapped in antinodes of a stand- 
ing optical wave. Each cloud is a 2D system. The number of 
atoms per spin state per square micron is coded in the tones 
of gray. The meaning of the tones is shown on the right. 



~ 0.1 Ep, where Ep is the Fermi energy. The interactions 
are tunable by means of the Fano-Feshbach resonance. 
This system can be used for qualitative and quantitative 
tests of many-body theories in two dimensions. 

In experiment, we use the atoms of lithium-6 in the two 
lowest-energy spin states, |1) and |2), with about equal 
populations: 50 ± 2.5%. In the limit of high magnetic 
field, these states have the same electronic spin projec- 
tion: rrij = —1/2 in the magnetic-field basis, and differ- 
ent projections of the nuclear spin, mi = 1 and mj = 0, 
respectively. When modeling condensed matter systems 
in experiments with cold atoms, states |1) and |2) may 
be regarded as analogs of the electronic spin- up and spin- 
down states, respectively. 

The standing-wave dipole trap is formed by two fo- 
cused, counterpropagating Gaussian beams with overlap- 
ping foci. The beams have identical power and polariza- 
tion, and the wavelength of 10.6 /im. The period of the 
dipole potential is 5.3 /im, which can be resolved on im- 
ages taken in the light at 671 nm wavelength, which is 
resonant to the lithium atoms. 

The degenerate Fermi gas is prepared following the 
procedure of Ref . [33j . During the first 6 s of preparation, 
10 8 — 10 9 atoms are collected in a magneto-optical trap 
(MOT) from an atomic beam. The standing- wave optical 
dipole trap is spatially overlapping with the MOT and is 
on during the MOT loading. After the MOT fields are 
turned off, ~ 10 6 atoms remain trapped in the standing- 
wave dipole trap whose depth is ~ 230 /iK. In order to 
start evaporative cooling, the s-wave scattering length 
is increased to a large negative value a — —3950 Bohr 
by switching on a nearly uniform magnetic field of 1020 
G in the — y direction. This value of the magnetic field 
corresponds to the Fermi side of a broad Fano-Feshbach 



resonance in s-wave scattering 3J]. During 1 s, the gas 
evaporates freely in a stationary potential. Afterwards, 
one of the beams forming the standing wave is gradu- 
ally turned off over 0.2 s. As a result, the gas adiabat- 
ically reloads to a cigar-shape dipole trap formed by a 
focus of a single traveling wave. During reloading the 
gas keeps evaporating and cooling. Further cooling is 
done by means of forced evaporation [35J during 10.6 s. 
Over this time, the trap depth is decreased by a factor of 
100 following the law U x {t) = (60 fiK) x [1 - t/(14 s)] 324 
by decreasing the trapping-beam power. In the single- 
beam dipole trap at low depth, the axial confinement is 
dominated by the small curvature of the magnetic field, 
which compresses the cloud in the axial direction. At 
this point, the second beam is reestablished during 0.2 s, 
which reloads the gas back into the standing-wave trap. 
Forced evaporation is then continued by decreasing the 
depth exponentially for 3.5 s by a factor of 10. After- 
wards, the trap depth is kept stationary for 0.5 s. Then 
the power of the trapping beams adiabatically, over 2.4 
s, increases to get a potential of desired height. This 
completes the preparation. 

In situ images of the gas are obtained by the absorp- 
tion imaging technique [16j : The atoms are irradiated by 
a 6 fj,s pulse of a uniform laser beam resonant to a cy- 
cling two-level transition for one of the two spin states. 
The imaging beam is shed in the y direction, which is 
opposite to the magnetic field and perpendicular to the 
axis of the trap cylindrical symmetry z. The shadow, 
which the atoms make in the imaging beam, is projected 
and recorded on a CCD camera. From the shadow, we 
reconstruct the column density distribution n,2(x,z) ex- 
actly accounting for saturation effects [33j due to finite 
intensity of the imaging beam: I = 2.6 mW/cm ~ ^sat- 
in counting the atom number, we correct for the effects 
that reduce the visible number of atoms: (i) Dopplcr shift 
caused by the acceleration due to the light-pressure force, 
(ii) diffractive spreading of the image of narrow clouds 
slightly beyond the objective size, (iii) small coupling of 
the imaging transition to a dark state, and (iv) fluores- 
cence. These effects reduce the apparent atom number 
by 18%, 8%, 3%, and 3% respectively. The spatial res- 
olution of our video system is 1.2 — 2.0 fjja.. The ob- 
served density distribution is shown in Fig. [21 The noise 
is dominated by the photon shot noise of the imaging 
beam. Each cloud in the figure represents an isolated 
two-dimensional Fermi system. 

The confining potential is a combination of a deep op- 
tical lattice in the z direction and Gaussian-shape con- 
finement in the transverse plane: 



V s (x) = sE r 1 — exp 



moj\{x 2 + y 2 ) 



2sE r 



cos 2 kz 



(2) 

where E r = h 2 k 2 /2m is the recoil energy [k — 
27r/(10.6 /jm)] and s is the dimensionless lattice depth. 



In formula ([2]) for the potential shape, we neglect beam 
divergence because the Rayleigh length zr ~ 5 mm 
is much bigger than the 400 /im long region where 
the atoms are prepared. At the bottom of each well, 
the potential is nearly harmonic as in Eq. ([T]). We 
use a version of the parametric resonance method 36 1 
to measure the frequencies uj±/2tt = 102 ±4 Hz and 
uj z /2it = 5570 ± 100 Hz and the dimensionless depth 
s = (hoj z /2E r ) 2 = 86.5 ± 3.0. The absolute trap depth 
value is sE r — (4.65 ± 0.18)fiw z = 1.23 fiK. Each trap is 
strongly anisotropic: uj z /uj± = 54.6 ± 2.4. 

Reduced dimensionality of the gas in a single cell can 
be proven by satisfying both of the following conditions: 
(i) The absolute majority of atoms is populating the ax- 
ial ground state and (ii) tunneling between the wells is 
negligible. 

The suppression of tunneling is assured by the lattice 
depth s = 86.5 » 1. The width of the lowest Bloch 
band is 1.1 x 10 _7 fiw 2 , which gives the tunneling time 
of ~ 260 s. For times much shorter than this, the gas 
remains kinematically two-dimensional. 

The population of the excited axial states may come 
from the Fermi statistics, thermal excitations, and mean- 
field interaction. In estimating the role of the statistics 
and temperature, we use the model of noninteracting 
Fermi gas in a parabolic potential. 

At zero temperature, due to the Pauli exclusion prin- 
ciple, atoms occupy energy levels up to the Fermi energy 
Ep = huj±^/2N, where N is the number of atoms per 
spin state in each cell and Ep is counted from the axial 
ground state. The necessary condition of dimensional- 
ity two is then Ep < fno z . We find N — 660 ± 60 by 
integrating the column density in each cell, where the 
error margins include deviations in the measurement be- 
tween different repetitions of the experiment as well as 
fluctuations between different clouds on the same pho- 
tograph (Fig. [5J. For the analysis we use 21 central 
clouds which are about equally populated. Therefore, 
E F = Huj ± V2N = 180±10nK= (0.67 ± 0.04) hw z : i. e., 
the Pauli exclusion principle does not create population 
of the axial excited state. 

The temperature is found by fitting the model one- 
dimensional density profile to the data. Figure [3] shows 
the experimental density distribution obtained by inte- 
grating n,2{x, z) along z in a single cell and averaging 
over 21 central cells. Also for the data of Fig. [3l trav- 
eling average over adjacent pixels is done to reduce the 
photon shot noise effect. The model density profile is the 
finite-temperature Thomas-Fermi distribution 
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where Li 3 / 2 is the polylogarithm function of order 3/2 
and \i is the chemical potential found from the condition 
N = J dxn\(x). The fitting parameters are T and wj_. 
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FIG. 3. One-dimensional density profile in a single cloud ob- 
tained by integrating 3D density distribution along y and z. 
The dots are the data averaged over the 21 central clouds. 
The black solid curve is the fit of formula © to the data. 
The blue dashed curve is the fit of a Gaussian distribution. 



We make tu± floating in the fit despite the availability of 
a measured value in order to account for the presence of a 
small mean field. The fitted curve is shown in Fig. [3] Fit- 
ting of a Gaussian (blue dashed curve), which would cor- 
respond to a nondegenerate gas, gives significantly larger 
deviation from the data. Fitting of the Thomas-Fermi 
profile © yields T = (0.10 ± 0m)E F = 18 nK. At this 
temperature, just 0.01% of the atoms are thermally ex- 
cited out of the axial ground state. 

The mean field can be treated as a perturbation be- 
cause the 3D interaction parameter kpa = —0.43 is small 
(kp = y / 2mEp/h) and because the system is far from the 
predicted geometric resonance [33|. We estimate the de- 
pletion of the axial ground state by using a simplified 
Hamiltonian for relative motion of two colliding particles 
in the center-of-mass reference frame: 



H 



Pi fhuj 2 z Z 2 
2^ + 2 



2irh 2 a (ri2(x, y)) 



8{Z). (4) 



where m = m/2 is the reduced mass; P z and Z are the 
operators of the relative momentum and coordinate of 
the two interacting atoms, respectively; 5 is the Dirac 
delta function; and (ri2(x,y)) = mEpf&ixY? is the trap- 
averaged 2D density distribution in the transverse plane. 
In this approximation, only the motion along z is ac- 
counted for. The last term accounts for the contact in- 
teraction, which depletes 0.2% of the atoms out of the 
axial ground state. 

As a result, the total depletion of the axial ground state 
is estimated to be ~ 0.2%. Together with the negligible 
tunnel rate, this proves that the axial motion is "frozen 
out," and we observe a series of isolated 2D Fermi sys- 
tems. 

The technique of trapping fermionic atoms in antin- 
odes of a standing wave has been used in the beautiful 



experiments devoted to molecular formation [38j, inter- 
ferometry [39|, and study of collisions [4fJ, |41|. None 
of these experimental systems have been found to be in 
the 2D Fermi-degenerate regime. In particular, in ex- 
periment |40j, the gas is not two-dimensional because 
the preparation process inevitably populates at least the 
two lowest Bloch bands and also because the tunnel 
time is shorter than the interval between preparation 
and probing. In experiment 4l|, the dimensionality of 
the lowest-energy samples is not clear and represents 
an interesting theoretical question. On one hand, the 
Fermi energy is above the first axial excited state. On 
the other hand, due to the strong interactions, the zero- 
temperature chemical potential could be still below Hu} z . 
The thermal populations of the excited axial states are 
unknown but can, in principle, be determined from the 
reported data via an appropriate theoretical analysis. To 
our knowledge, such analysis is not available at present. 
If in experiment 41| the thermal populations are neg- 
ligible, that system might be in the interesting quasi- 
two-dimensional regime, where the axial excited states 
are populated by the effect of the strong interactions 



alone [42 1 



In conclusion, we have prepared and directly observed 
a two-dimensional Fermi gas of atoms. This system may 
be a starting point for exploration of 2D Fermi physics 
and critical phenomena in a defect-free, controllable en- 
vironment. 
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